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Abstract. The kinetic energy distribution as a function of mass of final fragments 
(to) from low energy fission of '^^'^U, measured with the Lohengrin spectrometer by 
Belhafaf et ai, presents a peak around ra = 108 and another around to = 122. The 
authors attribute the first peak to the evaporation of a large number of neutrons 
around the corresponding mass number; and the second peak to the distribution of the 
primary fragment kinetic energy. Nevertheless, the theoretical calculations related to 
primary distribution made by Faust et al. do not result in a peak around m ~ 122. 
In order to clarify this apparent controversy, we have made a numerical experiment in 
which the masses and the kinetic energy of final fragments are calculated, assuming an 
initial distribution of the kinetic energy without peaks on the standard deviation as 
function of fragment mass. As a result we obtain a pronounced peak on the standard 
deviation of the kinetic energy distribution around m = 109, a depletion from to = 121 
to TO = 129, and an small peak around to = 122, which is not as big as the measured by 
Belhafaf et al. Our simulation also reproduces the experimental results on the yield of 
the final mass, the average number of emitted neutrons as a function of the provisional 
mass (calculated from the values of the final kinetic energy of the complementary frag- 
ments) and the average value of fragment kinetic energy as a function of the final mass. 
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1. Introduction 

One of the most studied quantities to understand the fission process is the fission 
fragment mass and kinetic energy distribution, which is very closely related to the 
topological features in the multi-dimensional potential energy surface of the fissioning 
system P . Structures on the distribution of mass and kinetic energy may be interpreted 
by shell effects on that potential energy, determined by the Strutinsky prescription and 
discussed by Dickmann et al. [2] and Wilkins et al. [3]. 

In order to investigate the dynamics of the fission process, the distribution of final 
fragment kinetic energy (e) as a function of final fragment mass (m), from thermal 
neutron induced fission of was measured by Belhafaf et al. [1], using the Lohengrin 
spectrometer. This distribution was represented by the mean value of kinetic energy 
(e) and the standard deviation (SD) (jg as function of m. The results present a first 
peak on cre(m) around m = 108 and a second one around m = 122, see Fig. [TJ 
The authors attribute the first peak to a large number of evaporated neutrons {v) 
around the corresponding primary mass (^4). Based on the small number of emitted 
neutron measured around A = 122, the second peak is attributed to the distribution 
of the primary fragment kinetic energy {E). However, theoretical calculations made by 
Faust et al. [5] do not result in a peak in SD for the distribution of primary fragment 
kinetic energy aE around A = 122. 

In order to clarify this apparent controversy, it is crucial to find the relation between 
the primary and the final kinetic energy distributions; the relation between the primary 
{Y{A)) and the final {Y{m)) mass yield; as well as the relation between the average 
value of the number of emitted neutron {u = z>) as a function of the primary fragment 
mass and the values corresponding to the experimental results. To address this question, 
we present Monte-Carlo simulation results for thermal neutron induced fission of ^33^/ 

1. e low energy fission of "^^^U 

2. Monte Carlo simulation model 

2.1. Fragment kinetic energy and neutron multiplicity 

In the process of thermal neutron induced fission of ^'^^[/, the excited composed nucleus 
"^^^U* is formed first. Then, this nucleus splits in two complementary primary fragments 
having Ai and A2 as mass numbers, and Ei and E2 as kinetic energies, respectively. 
Using relations based on momentum and energy conservation, the total kinetic energy 
of complementary fragments is. 




1 



(1) 



The total excitation energy is given by 



TXE = g + e„ - TKE, 



(2) 
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Figure 1. SD of the final fragment kinetic energy distribution as a function of the final 
mass m (• and *) as measured by Belhafaf et al. [4], and Faust et at [5], respectively; 
and SD as a function of primary mass (o) as calculated by Faust et al. [5] 



where Q is the difference between fissioning nucleus mass and the sum of two 
complementary fragments masses, and e„ is the separation neutron energy of "^^"^U. 
From Eq. ([T]) and ([2]), taking into account that ^41 + ^42 = 234, results 

234 

TXE = Q + e„ -E, (3) 

^ 234 - A ' ^ ^ 

where A and E are the primary mass number and kinetic energy, respectively, of one of 

the two complementary fragments. It is reasonable to assume that the excitation energy 

of one complementary fragment (E*) is proportional to the total excitation energy, thus, 

234 

E* oc TXE = g + e„ - ^^j^^^, (4) 

It is also reasonable to assume that the number (u) of neutrons emitted by a fragment 
is proportional to its excitation energy, i.e. 

u oc E*. (5) 

From relations (jl]) and we obtain a linear relation between u and E: 

u = a + bE,. (6) 

Taking into account that there is no neutron emission (z/ = 0) for fragments having the 
maximal kinetic energy [E = E^ax) and assuming that u = i? for the average value of 
fragment kinetic energy, the relation becomes: 

Emax E 

Let the parameter P define the maximal value of kinetic energy by the relation 

Emax = E + (8) 



(7) 
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Then, the relation ([7]) may be expressed as 

u = u{l~(3{^^)), (9) 

Because the neutron number N is an integer, it will be defined as the integer part of 
©, i.e. 



N = Integer part of (a + z/(l - P{ ))), (10) 

where a is used to compensate the effect of the change from a real number u to an 
integer number A^. 

2.2. Simulation process 

In our Monte Carlo simulation the input quantities are the primary fragment yield (Y), 
the average kinetic energy (E), the SD of the kinetic energy distribution ((Te) and the 
average number of emitted neutron (z/) as a function of primary fragment mass (A). 
The output of the simulation for the final fragments are the yield (Y), the SD of the 
kinetic energy distribution ((Te) and the average number of emitted neutron (z/) as a 
function of final fragment mass m. 

For the first simulation, we take Y and E from Ref. The first SD curve 
is an extrapolation of calculation results obtained by Faust et al. p]. Then, we adjust 
Y{A), vi^A), E{A) and aE^A) in order to get y(m), i/, e(m), aeijn) in agreement to 
experimental data. 

In the simulation, for each primary mass A, the kinetic energy of the fission 
fragments is chosen randomly from a Gaussian distribution 

{E-E) 



P{E) = — — exp 
\/1t{Oe 



2a% 



11) 



where P{E) is the probability density of energy with mean value E and SD aE- 

For each E value, the simulated number of neutrons N is calculated with the relation 
(ITOl) . The final mass of the fragment will be, m = A — N. Furthermore, assuming that 
the fragments loose kinetic energy only by neutron evaporation and not by gamma 
emission or any other process, and neglecting the recoil effect due to neutron emission, 
the kinetic energy e(m) of the final fragment will be given by 

e(m) = (1 - ^)E, (12) 

With the ensemble of values corresponding to m, e and A^, we calculate Y{m), e(m), 
aeijn) and vim). 

To obtain an acceptable statistics during the simulation, we have considered a total 
number of fission events of ^^^[/ of the order of 10*, and we have computed the SD of 
all the relevant quantities by means of the following expression: 
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Figure 2. Simulation results for the primary (A) and final (o) mass yields are 
presented together with experimental data (•), taken from Ref. [3] 
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Figure 3. The average number of emitted neutrons from fission of ^'^'^C/ : as a function 
of the primary fragment mass A (A), as a function of final fragment mass (©) both as 
result of simulation and experimental data (•), taken from Ref. [6j 

where e(m) is the mean value of the kinetic energy of final fragments with a given mass 
m, and Nj{m) is the number of fission events corresponding to that mass. 

3. Results and interpretation 

The simulated final Y{m) and the primary Y{A) mass yield curves are illustrated in 
Fig. O As expected, due to neutron emission, the Y{m) curve is shifted from Y{A) 
towards smaller fragment masses. The simulated average number of emitted neutron 
z>(m) curve is shifted from iy{A) in a similar way as Y{m) relative to Y{A){see Fig. [3]). 

As stated in sect. |2], the primary kinetic energy {E{A)) is generated from a Gaussian 
distribution, while the final kinetic energy (e(m)) is calculated through Eq. (fT2l) . 

The plots of the simulated mean kinetic energy for the primary and final fragments 
as function of their corresponding masses are shown in Fig. |H In general, the simulated 
average final kinetic energy curve as a function of final mass (e(m)) display a shift 
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Figure 4. Mean kinetic energy of the primary fragments (A) and mean kinetic energy 
of the final fragment (©), as a resuh of simulation in this work, to be compared to 
experimental data (•) as measured by Belhafaf et al. f4j. 
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Figure 5. Assumed SD of primary fragment kinetic energy distribution (A) and 
simulated in this work final fragment kinetic energy distribution (©), to be compared 
to results of calculations (o) made by Faust et al. [5^ and experimental data (•) as 
measured by Belhafaf et al. [4j 



roughly similar to that of Y{m) curve, and a diminishing of e relative to E values given 
by relation (fT2|) with N = u. 

Furthermore, Fig. [5] displays the SD of the kinetic energy distribution of the primary 
fragments {aE{A)) and the SD of the kinetic energy of the final fragments ((Je(m)). The 
simulated initial distribution of the kinetic energy have no peaks on the SD. The plots 
of (Te(m) reveal the presence of a pronounced peak around m = 109 in agreement with 
the experimental results obtained by Belhafaf et al. [1] and Faust et al. [5] , respectively. 
The peak on the SD around m = 122 resulting from simulation is not as big as the 
obtained by Belhafaf et al. Moreover a depletion on the SD in the mass region from 
m = 121 to m = 129 is obtained as a result of simulation. The simulated primary 
fragment kinetic energy distribution (see Fig. OA) does not present peaks in the range 
of fragment masses A from 90 to 145. If one simulates an additional source of energy 
dispersion in erg, without any peak, no peak will be observed on Ce. Both the shape and 
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Figure 6. SD of final fragments kinetic energy distribution calculated assuming that 
i) Y(A) is constant, ii) (Je{A) is constant, iii) fragments with E > E do not emite 
neutrons and fragments with E < E emit one neutron and iv) neutron emission have 
no recoil effect on fragment kinetic energy. E{A) values are taken from data. 



height of the peaks of ae{m) are sensitive to the value of parameter a and /9 appearing 
in Eq. (JTOl) . A higher value of a will produce a larger peak of SD. The effect of on peak 
depends much on mass region. For the region m = 108, a higher value of (3 will produce 
a larger peak of CTg. The simulated results for ae{m) presented in Fig. [5] were obtained 
with a = 0.62 and /3=0.35. The presence of a peak at m = 108 could be associated 
with neutron emission characteristics (approximately i? = 2) and a very sharp fall in 
kinetic energy from E =96 MeV to E =90 MeV, corresponding to ^4=109 and ^=111, 
respectively. A similar result was obtained for low energy fission of ^se^/ 

In order to more easily interpret the influence of variation of E{A) on cre(m), we 
derive an analytical relation assuming that: i) Y(A) is constant, ii) (Te{A) is constant 
and iii) fragments with E > E do not emit neutrons and fragments with E < E emit one 
neutron. Then for each final mass there is a contribution from fragments with primary 
mass m that do not emit any neutron and from fragments with primary mass m + 1 
that emit one neutron. With these conditions we can show that. 




(14) 



where AE = E{m + 1) — E{m) 

The equation (HM corresponds to a parabola with a minimum value 



which occurs when 

AE = 2J-aE = l.Q(TE, (16) 

V 71 

As we can see on Fig. [6l the aeijn) curve, calculated with relation (IT^ . presents 
a peak around m = 109 in reasonable agreement with the experimental data. In that 
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Figure 7. SD of final fragments kinetic energy distribution calculated under the 
assumption thar i) (T£;(A) — 5 MeV , iii) E{m + 1) = E{m) and iv) neutron emission 
have no recoil effect on fragment kinetic energy. Y{A) values are taken from data. 



region Ai? < 0, then from relation (|T^ it follows that o"e(m) > a"s(y4). The depletion 
on the simulated a"e(m) on the mass region between m = 121 and m = 129 is explained 
by the fact that AE > 0. Using relation (fT^ . we obtain that aeijn) < aE^A). 

In order to more easily evaluate the influence of the variation of Y on aeim), we 
derive an analytical relation assuming that ( i) Y{m + 1) = rY{m), (ii) aE{A) are 
constant, (ii) E{m + 1) = E{m) and (iii) neutron emission have no recoil effect on 
fragment kinetic energy. Then we can show that 



(Te(m) = Ge 



1-1 



1 



(17) 



The ae{m) values calculated with relation (|T71) are lower than aE{A) and higher 



than 



ae{m) 



TT 



-aE = 0.6(Te, 



The SD curve calculated with this relation is presented in Fig. [71 Using relation 
(fTSll we can estimate a peak at m = 122 assuming that around this mass Y increases 
very rapidly with A except that y(123) = y(122). However, we can not reproduce the 
pronounced peak obtained by Belhafaf et al. J4j. 



4. Conclusions 



Using a simple model for the neutron emission by fragments, we have carried out a Monte 
Carlo simulation for the mass and kinetic energy distributions of final fragments from 
thermal neutron induced fission of "^^^U. In comparison with the primary fragments, 
the final fission fragments have eroded kinetic energy and mass values which give rise 
to the appearance of peaks in the SD of the final fragments kinetic energy as a function 
of mass cre(m) around m = 108. This peak, which agrees with experimental results 
obtained by Belhafaf et al. |4j and Faust et al. [5] is a consequence of neutron emission 
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and variations of primary fragments yield {Y(A)) and mean kinetic energy E[A) curves. 
As a result of the simulations one obtains a depletion on the SD in the mass region from 
m = 121 to m = 129, and an small peak on the SD around m = 122 which is not as 
big as the measured by Belhafaf et al. Our simulation also reproduces the experimental 
results on the yield of the final mass. 
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